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The electromagnetic decays of the S°(1385) and A(1520) hyperons were studied in photon-induced 
reactions 7 p ^ A''^A(1116)7 in the CLAS detector at the Thomas Jefferson National Accelerator 
Facility. We report the first observation of the radiative decay of the E°(1385) and a measurement 
of the A(1520) radiative decay width. For the E°(1385) ^ A(1116)7 transition, we measured a 
partial width of 479 ± 120(stat)l®QQ(sys) keV, larger than all of the existing model predictions. For 
the A(1520) ^ A(1116)7 transition, we obtained a partial width of 167 ± 43(stat)^^2(sys) keV. 

PACS numbers: 14.20.Jn,13.30.Ce,13.40.Hq 


INTRODUCTION 

The low-lying neutral excited-state hyperons E°(1385), 
A(1405), and A(1520) were discovered in the 1960s, but 
their quark wave functions are still not well-understood 
and experimental studies of their properties have been 
scarce since the early 1980s. The electromagnetic decays 
of baryons produced in photon reactions provide an es¬ 
pecially clean method of probing their wave functions. 
Baryons with a strange quark have an additional degree 
of freedom which aids in the study of multiplet mixing 
and non 3-quark admixtures. Recently there has been a 
renewal of interest in this held, e.g. electro-production of 
the A(1520)Q. This paper reports the results of a non¬ 
model dependent measurement of the radiative decay of 
the E°(1385) and A(1520). 

The non-relativistic quark model (NRQM) of Isgur and 
Karl|3 has been remarkably successful in predicting the 
masses and widths of N* and A* states, but less suc¬ 
cessful in the strange sector. Several competing models 
for hyperon wave functions have been proposed. Mea¬ 
suring the transitions Y A(1116)7 and Y —s- E(1193)7 
provides a means of differentiating between these mod¬ 
els. Calculations have been done in the framework 
of NRQM0, a relativized constituent quark model 
(RCQM)|3, a chiral constituent quark model (yCQM) 
that includes electromagnetic exchanre currents between 
the quarks 0, the MIT bag modelQ- the chiral bag 
model| 3 , the bound-state soliton model@|, a three-flavor 
generalization of the Skyrme model that uses the collec¬ 


tive approach instead of the bound-state aDDroachPll(i|. 
an algebraic model of hadron structure|H| , heavy baryon 
chiral perturbation theory (HByPT)| 12|. and the l/N^ 
expansion of QCD[iH. The radiative widths in keV are 
tabulated in Table ^ The A ^ pj width is included for 
comparison. 

The photon decay spectrum of the low-lying excited 
state hyperons is shown in Fig. ^ The widths given in 
Table m can be qualitatively estimated using SU(3) sym¬ 
metry. The A(1116) and the E°(1193) are in the S=l/2+ 
SU(3) octet and the E°(1385) is in the S=3/2+ SU(3) 
decuplet. The A(1116) has the two light quarks in the 
s orbital in a spin S=0, isospin T=0 configuration. The 
E°(1193) and the E°(1385) have the light quarks in a 
spin S=l, T=1 configuration. All three hyperons have 
the strange quark in the s orbital. Decuplet to octet ra¬ 
diative decays are dominated by an Ml transition with 
a spin-flip of one quark. The SU(3) model prediction 
of the ratio E* ^ E 7 to the E* ^ A 7 is ~ g times 
kinematic factors. This, plus the fact that most of the 
constituent quark model calculations 13.1^ IR listed 
in Table Q] used the impulse approximation, leads to a 
very narrow range of predictions, (265-273 keV) for the 
E* —> A 7 reaction and (17.4-23 keV) for the E* ^ E 7 
reaction. The A(1405) and A(1520) have light quarks 
in the s orbital with S=0, T=0 and the strange quark 
in a pi and p| orbital respectively. The radiative de¬ 
cays A(1520) ^ E*, E and A(1405) ^ E*, E require that 
the strange quark make a transition from a p orbital to 
an s orbital with a simultaneous spin flip of one of the 
light quarks. These transitions are thus forbidden by the 
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FIG. 1: Photon decay spectrum of low lying excited state hy- 
perons. The transitions shown as dashed lines are suppressed. 


one body nature of the electromagnetic operator. They 
can proceed only via configuration mixing introduced by, 
e.g. the QCD hyperfine interaction, which leads to a 
wider range of model predictions. This is explained in 
more detail in an excellent review of the experimental 
and theoretical situation in M 

Experimental measurements have been sparse. The 
results are tabulated in Table Q] The Afl520) Ay 
transition has been measured by Mast et nl. jl (t|| using 
a K~ beam with a liquid-hydrogen bubble chamber, by 
Bertini et al. (unpublished) with a liquid-hydrogen 
target viewed by a Nal detector, and by Antipov et aL|l^ 
using a high-energy proton beam on carbon and copper 
targets. Antipov et al. measured the K~^, p and tt” in a 
magnetic spectrometer and detected the decay photons 
using an electromagnetic calorimeter. These are the only 
direct measurements in the literature. Burkhardt and 
Lowe[l^ extracted model dependent branching ratios for 
A(1405) radiative decay from the kaon-proton capture 
data of Whitehouse et al. The radiative decay of the 
E°(1385) has never been observed (Meisnerj^ reports 
one event); only upper limits for the branching ratios 
have been established|2^. 


EXPERIMENT 

In the current experiment, the low-lying excited-state 
hyperons were studied in the reaction ■jp —> K^p'K~ X us¬ 
ing the CEBAF Large Acceptance Spectrometer (CLAS) 
in Hall B at the Thomas Jefferson National Accelera¬ 
tor Facility. The data were from the GIC running pe¬ 
riod September to October 1999. The primary elec¬ 
tron beam was converted to a photon beam with a 
thin radiator of 10“^ radiation lengths. The scattered 
electron is momentum-analyzed by a photon tagging 
spectrometer [ 2 ^ with a resolution of AE/E = 10“^. Pho¬ 
tons were tagged over a range of 20-95% of the inci¬ 
dent electron beam energy. The electron beam energies 
were 2.445 GeV, 2.897 GeV, and 3.115 GeV and the cur¬ 
rents were typically 6 nA. The target was liquid hydrogen 
in a cylindrical cell of 17.9 cm length and 2 cm radius. 
The GLAS detector consists of six individually instru¬ 
mented segments, each consisting of three layers of drift 
chambers and a shell of 48 time-of-flight scintillators. Six 
superconducting magnets provided a toroidal magnetic 
field, with negative particles bent toward the beam di¬ 
rection. The trigger consisted of a triple coincidence be¬ 
tween the photon tagger, the time-of-flight system, and 
a small scintillation detector (the “Start Counter”]^) 
surrounding the target scattering chamber. Only one 
charged particle in the GLAS was required in the trigger 
to accommodate the 6 experiments that were running si¬ 
multaneously. A total of 1420M triggers were collected 
at 2.445 GeV, 845M at 2.897 and 2280M at 3.115 GeV. 

Particle Identification 

Gharged hadrons were identified using momentum and 
time-of-flight information. The processed data files were 
filtered for events containing one K'^, one 7r“, and one 
proton track in coincidence with the incident tagged pho¬ 
ton. Kaon candidates were chosen using a broad range in 
mass (0.35-0.65 GeV). The 7r“ candidates were selected 
with a mass range of <0.3 GeV and proton candidates 
with a range of 0.8-1.2 GeV. A minimum momentum cut 
of 0.3 GeV/c was applied for the kaons and protons and 
0.1 GeV/c for the pions. The hadron mass spectrum for 
events that survive the filter is shown in Fig. The 
kaon peak sits on top of a large background due to high 
momentum pions with poorly determined mass. 

To further refine the kaon identification the difference 
At between the time at the target for the kaon candidate 
and what it would be for a true kaon was computed: 

A( = {troF - U„) (^1 - 

AC, =-A 


( 2 ) 
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A(1232) 

E°(1385) 

A(1405) 

A(1520) 

Model 

P7 

A(1116)7 

E"(1193)7 

A(1116)7 E'’(1193)7 

A(1116)7 

E"(1193)7 

NRQM[;i, J 

360 []J 

273 

22 

200 

72 

156 

55 

RCQM[^ 


267 

23 

118 

46 

215 

293 

xCQM[a 

350 

265 

17.4 





MIT BaE[3l 


152 

15 

60, 17 

18, 2.7 

46 

17 

Chiral Bag[^ 




75 

1.9 

32 

51 

Soliton[^ 


243, 170 

19, 11 

44, 40 

13, 17 



Skyrme[9,J^ 

309-348 

157-209 

7.7-16 





Algebraic modeUll] 

343.7 

221.3 

33.9 

116.9 

155.7 

85.1 

180.4 

HBxPTfljT 

(670-790) 

290-470 

1.4-36 



1/Nc expansion[^ 


298 ± 25 

24.9 ±4.1 



Previous Experiments 

640-720130] 

<2000 [22.1 

<1750[22] 

27±8|49j 

10±4[19j 

33±11[17] 

47±17fl7] 






23±7[12j 

134±23ljJ 








159±33±26|T2j 


This experiment 


479 ± 120+100 


167 ± 43+12 


TABLE I: Theoretical predictions and experimental values for the radiative widths (in keV) for the transitions Y A(1116)7 
and Y —> E( 1193)7. Some models have mnltiple predictions that depend on different assumptions. For comparison the 
predictions and experimental value are quoted for the A(1232) —> p 7 transition. ^The results for HBxPT[T 3 are normalized 
to the quoted empirical range (in parentheses) for the A ^ pj transition. 



FIG. 2: Particle Identification: Hadron mass from TOF and 
momentum information multiplied by the sign of the charge 
of the particle. The shaded curve is the mass spectrum after 
the PID cuts. 


where troF — tvert is the flight time between the nterac- 
tion vertex and the he Time-of-Flight array. We require 
|Af| ^ 0.G7 ns. Since the experiment consists of tv^o 
physically separate systems, the Tagger and the CLAS 
detector, we require that the time at the interaction ver¬ 
tex measured by the two systems agree to within 5% of 
the flight time between the Start Counter and the TOF 
paddles. 

The CLAS detector does not cover the full angular 
range in 6 or (j). Some angular regions are shadowed by 
the toroidal coils. The shadow region broadens in (j) as 
a function of decreasing 0 as seen from the center of the 
target. All tracks were required to be in the region of 
well-understood acceptance by applying a fiducial cut of 


the form 


6 » > 4.0-b 


510.58 

(30 - 0)1-5518 ’ 


(3) 


where 4> is the azimuthal angle folded onto the range 0- 
30°. We also require (j) < 26°. 

Some of the ”kaon” events are really misidentified 7r+. 
This can be seen in Fig. 13 where all events are plotted 
assuming that all kaon candidates are really misidentified 
TT'*" and compute the missing mass squared for the reac¬ 
tion 7 p ^ p7r'''7r“(A). The prominent spike at zero mass 
squared indicates jp contamination and a 7r° 

peak is clearly evident but at a much reduced level. The 
expected distribution for good K'^ events goes to zero 
for zero p'k~^tt~ missing mass squared. We require the 
missing mass squared from this calculation to be greater 
than 0.01 GeV^ to eliminate for example p —> tt+tt" con¬ 
tamination. We did not cut above the 7r° peak in Fig. |3 
because that would have cut into the good K'^ events. 
The hadron mass spectrum after all of the above cuts 
have been applied are shown as the shadowed histogram 
in Fig. [3 

The kaon momentum is corrected for average dE/dx 
losses in the target material, target wall, carbon epoxy 
pipe and the start counter depending on the position of 
the primary vertex, which is approximated by the inter¬ 
section of the proton and kaon tracks. The ground-state 
A is sufficiently long-lived that it decays a measurable 
distance from the primary vertex. The secondary ver¬ 
tex is determined by the intersection of the proton and 
7 r“ tracks. The proton and tt~ tracks are corrected for 
average dE/dx losses according to the position of the sec¬ 
ondary vertex. 

The four-momentum of the A(1116) was reconstructed 























imize 


A) = (770 - rfj'^V ^{rfo - rf) + 2X^ f (4) 



FIG. 3: Pion contamination: Mass squared {Mx) for the 
7 P —> pn^TT~ {X) reaction where the was a potentially 
mis-identihed kaon. 



FIG. 4: A identihcation: proton-7r invariant mass 


from the proton and 7r“ four-momenta (Fig. 0. The 
Gaussian resolution of the A peak is about a = 1.3 MeV, 
consistent with the instrumental resolution. The excited- 
state hyperon mass spectrum for the region between 1.25 
GeV and 1.75 GeV requiring the pTr~ invariant mass to 
be in the range 1.112-1.119 GeV is shown in Fig. UK- Fig. 
Ep shows the mass Mx from the reaction yp ^ A(V). A 
clear peak at the mass of the AT* (892) is seen. The peak 
at the mass is due to accidentals under the TOF 
peak. This background is eliminated by requiring Mx > 
0.55 GeV. Fig. 0 shows the missing mass squared for 
the reaction yp —> Ar+A(A) after the foregoing cuts have 
been applied. A prominent peak shows up at and a 
smaller peak at zero missing mass squared. The counts 
above the 7r° peak are typically due to yp ^ K'^YP{X). 


Kinematic fitting 

A better approximation to the primary and secondary 
vertices can be found using kinematic fitting. We used 
the Lagrange multiplier method . The unknowns are 
divided into a set of measured variables {if) and a set of 
unmeasured variables {f) such as the missing momentum 
or the 4-vector for a decay particle. For each constraint 
equation a Lagrange multiplier \i is introduced. We min- 


by differentiating with respect to all the variables, lin¬ 
earizing the constraint equations and iterating. Here pb 
is a vector containing the initial guesses for the measured 
quantities and V is the covariance matrix comprising the 
estimated errors on the measured quantities. We iterate 
until the difference in magnitude between the current 
and the previous value is <0.001. The covariance matrix 
V for each track returned by the tracking code does not 
contain the effects of multiple scattering and energy loss 
in the target cell, the carbon epoxy pipe, or the start 
counter. To correct for this we apply multiple scatter¬ 
ing and energy loss corrections to the diagonal matrix 
elements. 

The first step in the fitting procedure is to fit the pro¬ 
ton and 7r~ tracks with the A hypothesis. This is a 2G fit. 
There are six unknowns (pA, Pv 2 ) and eight constraint 
equations. 


/ = 


Up + Att — Aa 
Pp+P-K - PA 

(y - y^)p% -{z- z^)py 

{x-x^)pI- {z- z^)pI 

{y - yp)pf -{z- zp)pl 

(x - Xp)pf -{z- Zp)pf, 


(5) 


The distribution for this fit is shown in Fig. [7|A and 
the Gonfidence Level plot is shown in Fig. [7j3. The curve 
is the result of a fit to the histogram using the function 
form of a x^ distribution with two degrees of freedom 
plus a flat background term. Explicitly, 


fix^) 


^p-P2X^l2 

2 


+ ^^3- 


( 6 ) 


The fit result {P 2 = 0.558) suggests that we are under¬ 
estimating the errors in the proton and tt~ tracks, but 
the shape is close to the expected shape. The Confidence 



FIG. 5: (A) Missing mass for the reaction yp —> K'^{X). (B) 
Missing mass for the reaction yp ^ A(A). 
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Level is given by the equation 


CL = 


f{z;n)dz 


(7) 


where f(z;n) is the probability density function with 
n degrees of freedom. 

The second step is to use these Kaon and Lambda 
tracks to obtain a better primary vertex. This is a 1C 
fit. There are 3 unknowns (iVi) and four constraint equa¬ 
tions. The distribution for this fit is shown in Fig. EP 
and the Confidence Level plot is shown in Fig. [7|I). The 
curve in Fig. w is the result of a fit to the histogram 
using the functional form of a distribution with one 
degree of freedom plus a flat background term. Explic¬ 
itly, 


/(X^) = 


Pi 


,-P2xV2 


V2F(1/2) 


Ps. 


( 8 ) 


with a fit result of P 2 = 0.507. We require the proba¬ 
bility of the A ^ pTT~ fit and the primary vertex fit be 
< 0.5% of exceeding for an ideal distribution. The 
improved kaon and lambda four vectors are used to com¬ 
pute the excited-state hyperon mass spectrum and the 
missing mass squared. 

Fig. compares the z-position of the primary vertex 
from the improved fitting procedure to the naive kaon- 
proton result. We apply a target z-position cut for the 
primary vertex between -10.0 cm and -1-9.0 cm and a 
radial cut of 2 cm. These cuts were chosen to ensure 
that the primary event came from the target region. The 
proper time of the A decay is plotted in Fig. 03. An 
exponential fit to the data gives a decay constant of 
7.62 ± 0.09 cm which is comparable to the PDG value 
of 7.89±0.06 cm. To verify that the target walls do not 
make a significant contribution to our yields, we applied 
the analysis procedure described above to the empty- 
target data. For the empty target runs the beam current 
ranged between 10 and 24 nA and averaged about 15 
nA. The results from analyzing about 33 million empty- 
target events (corresponding to approximately ^ of the 



FIG. 7: and Confidence Level distributions for the A ^ 

P'K~ fit (A and B) and the K vertex fit (C and D). 



FIG. 8: A. Z-position of primary vertex. Solid histogram: 
KA fit. Dashed histogram: Kp fit. B. Lambda decay proper 
time in units of cm. The excited-state hyperon mass was 
greater than 1.25 GeV for both plots. 



FIG. 6: Missing mass squared for the reactions yp —> 
K+k{X). 


target full integrated photon flux) are shown in Fig. O 
We obtained 25 A(1116) candidates within the proton-7r“ 
invariant mass range of 1.112-1.119 GeV (Fig. |5J4). The 
z distribution is shown in Fig. The hyperon mass 
distribution for those events satisfying the vertex cut is 
shown in Fig. EP. Figure EP shows the missing mass 
squared distribution for hyperon masses in the 1.34-1.43 
GeV range. There are no counts near zero missing mass 
squared and only two near rn^o- Both of these counts 
have z and r positions within the target volume. They 
correspond to interactions with the residual (cold) hy¬ 
drogen gas in the target. From this we conclude that 
the background due to interactions with the walls of the 
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FIG. 9: Empty target results. A. Proton piminus invariant 
mass. B. Vertex z-position. C. Hyperon mass. D. Missing 
mass squared for Mv=1.34-1.43 GeV. 

target cell is negligible. 

To achieve 7 / 7 r° separation, the events were sorted ac¬ 
cording to topology using kinematic fits with two hy¬ 
potheses 

Rl: 7 p —> itr+ATT^ 1C 
R2: 7 p^R:+A 7 1C 

The corresponding constraint equations are 

E_ ^beam 4“ Afp 74^ Ex _ Q /g\ 

[ Pbeam - PK - PA - PX \ ' ^ 

Here X is a missing 7 r° or a missing 7 . 

The distributions for reactions Rl and R2 are shown 
in Fig. HHK and nnp, respectively. The hyperon mass 
range was 1.25-1.75 GeV. The corresponding Confidence 
Levels plots are shown in Fig. IlOb and IIOL). For Rl 
we obtain the expected shape for a distribution with 
one degree of freedom. For R2 the x^ values indicates 
that the radiative decay hypothesis is inconsistent with 
most of the events. The dashed curve in Figure iint> is 
the Confidence Level for hypothesis R2 for those events 
which do not satisfy hypothesis Rl at the 5% level. We 
now see a shape consistent with a x^ distribution with 
one degree of freedom. 

Fig. HI] shows the missing mass squared distributions 
for a representative set of x^ cuts. For the purpose of the 
plot, we require xjn > 3.841 and x‘r 2 < 3.841 to isolate 
the radiative channel (case B). To isolate the pion chan¬ 
nel (case A), we require Xri < 3.841 and Xm — 3.841. 


FIG. 10: 8-nd Gonfidence Level distributions for the two 

reactions Rl (A and B) and R2 (C and D). The dashed curve 
in D is the R2 Gonfidence Level with the Rl reaction vetoed 
with x^=3.841. 

Case C is the “ambiguous” case where both xjn < 3.841 
and Xm < 3.841. Case D consists of those events that do 
not agree with either the radiative channel or the pion 
channel, for which Xm ^ 3.841 and xji 2 ^ 3.841. For 
a 1C fit = 3.841 corresponds to a 5% probability of 
exceeding x^ for an ideal x^ distribution. The “ambigu¬ 
ous” events are most likely to be XP events. 

Case D events are most likely be XP events. 

Fig. 1121 shows the corresponding hyperon mass spectra. 
Fig. II 21 \ is dominated by the E°(1385) ^ Att® chan¬ 
nel, for which the branching ratio is ^ 88 % [s^. We cal¬ 
culated the E(1385) radiative transition relative to this 
channel. The A(1520) peak shows up in Fig. 11 214 be¬ 
cause of the decay channels A(1520) —> E° 7 r° (BR=14%) 
and A(1520) ^ Atttt (BR=10%). 


Double Bremsstrahlung 

The 7 channel does not show the structure expected 
from hyperon photon decays. The structure was found 
to be masked by a background resulting from double 
bremsstahlung in the radiator. The reaction 71 -|- X 2 P 
K'^A + 7 i can mimic the reaction XP K'^ Ax ■ But 

in this case the missing momentum from the reaction 
XP K'^A{X) points along the +z direction (along the 
beam). This can also happen if the event is an accidental 
or inefficiencies in the tagger plane allow the wrong elec- 





























































FIG. 11: Missing mass squared distributions for x^iigh = FIG. 12: Hyperon mass distributions for x^high = Xlow = 
X%ow =3.841. Gases A-D are explained in the text. 3.841. Cases A-D are explained in the text. 


Iron to be selected. This problem is illustrated in Fig. El 
Fig. EK shows the off-z-axis momentum p]_ = 
for the candidate missing particle. 

This misidentification should happen for ground-state 
i;°(1193) production as well. A subset of the data filtered 
on the hyperon mass region between 1.0 and 1.25 GeV, 
was used to isolate E°(1193) events. The tr from a Gaus¬ 
sian fit to the peak is about 6.6 MeV, corresponding 
to a full width at half maximum of F = 2.354tT = 15.6 
MeV. This is a measure of the hyperon mass resolution. 
Apart from the hyperon mass range, the same set of cuts 
was used to analyze these data as for the excited-state 
sample. Fig. I13C shows the distribution in p\ for this 
data set. Figures 113b and EP compare the effect of 
two choices for the p\ cut on the hyperon mass distri¬ 
bution for the case where the 7 channel is favored. The 
histograms show the distributions in hyperon mass for 
those events that were cut out. Histograms lOtB andEP 
both look like exponentially falling distributions. Fig. El 
shows the corresponding hyperon mass spectra after ap¬ 
plying the p\ = 0.015 GeV^ cut. The histogram now 
shows the expected structure for the E(1385) ^ Ay and 
A(1520) ^ Ay reactions. Gomparison with I12IA shows 
that this cut also reduces the number of Att^ events seen. 
The Monte Garlo simulation □ is used to correct for this 
reduction. The p\ = 0.015 GeV^ cut will be used for the 
rest of the analysis. 



FIG. 13: Effect of cut on y channel. A. Thepx momentum 
spectrum for K'^A{X) events. B. The y channel cut distibu- 
tion for a 0.0004 GeV^ cut. G. The p\ momemtum spectrum 
for K'^T,°{X) events. D. The y channel cut distibution for a 
0.015 GeV^ cut. 
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FIG. 14: Hyperon mass distributions for Xlow = Xhigh = 
3.841 with px = 0.015 GeV^ cut. The labels are explained in 
the text. The yield of Att® and A 7 events in A and B were 
extracted by fitting the data with a relativistic Breit-Wigner 
(solid line) and a polynomial background (dashed line). In 
D the dashed histogram shows the contribution due to the 
E(1385) alone. The dotted histogram is the A(1405) con¬ 
tribution alone using the M-matrix parameterization for the 
shape. 


ACCEPTANCE 

A detailed Monte Carlo simulation of the CLAS de¬ 
tector was performed using GEANT 3.21 for each of the 
three electron beam energies. Table lTm lists the set of re¬ 
actions for which we generated events. The experimental 
photon energy distribution was used to determine the 
energies of the incident photons in the simulation. Rel¬ 
ativistic Breit-Wigner shapes were used for the E(1385), 
A(1520) and K* mass distributions. For the A(1520) the 
exponential slope for the t-dependence was 2.0 GeV~^. 
The angular distribution for the radiative decay of the 
A(1520) in its rest frame was taken to be proportional to 
5 — 3 cos^ 0 according to the result obtained by Mast, et 
al. . The same distribution was used for the Ett chan¬ 
nels and for the E(1385) decays. The model of Nacher, 
et al. [i^l with a flat angular distribution was used for 
the A(1405) decay channels. 

The incident photon energy dependence and t- 
dependence were adjusted to fit the data for the S(1385) 
reactions independently for each of the electron beam en¬ 
ergies. The data and MC were cut on the Y* mass range 
of 1.34-1.43 GeV and on the 7r° peak found in Fig. I14t \ 


FIG. 15: Momentum and angular distributions for MC 
(dashed histograms) and data (points with error bars) for the 
1.34-1.43 GeV hyperon mass region. 

to isolate the E(1385) ^ Att® channel. We plotted the 
ratio of the data/MC versus photon energy Ebeam- The 
resultant curve was fitted with a function of the form 
A/Ebeam + B/E^^^j^. We used this to modify the photon 
energy dependence of the E(1385) production cross sec¬ 
tion in the MG. The above procedure was then iterated. 
The exponential slope parameter was varied until the MG 
and data t distributions matched reasonably well. The 
exponential slope for the modified t-dependence was 1.0 
GeV~^. To check the quality of the simulation, we com¬ 
pared the momentum distributions for the Monte Garlo 
and the data for the kaon, proton, and pion tracks. The 
simulated events were analyzed with the same cuts de¬ 
scribed above. The results for the second iteration for 
the MG simulation are shown in Fig. ESI The agreement 
between the MG and the data for the pion, proton, and 
kaon momenta and the kaon lab angle is good. 

Fig. II til compares the data for the 1.49-1.55 GeV mass 
range and the missing mass squared in the range 0.018- 
0.075 GeV^ to the A(1520) ^ S°7r° Monte Garlo results. 
The MG results have been scaled by 0.185. The agree¬ 
ment between the MG and the measured momenta distri¬ 
butions is very good and the kaon angular distributions 
agree reasonably well. 

In order to check that the p]_ cut did not introduce a 
bias of the Monte Garlo results with respect to the data, 
we studied the yield of S(1385) ^ Att^ events in the 
data and the corresponding Monte Garlo. For the data 
we used the standard cuts and performed the same 
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FIG. 16: Momentum and angular distributions for MC 
(dashed histograms) and data (points with error bars) for the 
1.49-1.55 GeV hyperon mass region. 


pj cut (Gel/^) N(data) N(MG) N(data)/N(MC) 


0.005 

4021 

11037 

0.364T0.007 

0.010 

3500 

9860 

0.355T0.007 

0.015 

2878 

8148 

0.353±0.008 

0.020 

2191 

6191 

0.354±0.009 


TABLE II: Comparison of yields between the E(1385) —> Att® 
simulation and the data as a function of the p\ cut. The errors 
are statistical only. 


kind of fit to the hyperon mass distributions as described 
earlier. The hyperon mass range was 1.34-1.43 GeV. The 
results are tabulated in table HD The data and the MC 
yields agree as a function of the p\ cut. 

Table ITIII lists the acceptances for the case where 
X^iiGH — Xlow — 3.841. In the table At^ and refer 
to the fraction of surviving events relative to the number 
of thrown events that satisfy the Att^ and Ay hypotheses, 
respectively, and A-^.,^ refers to those events that do not 
satisfy either hypothesis. 


E°(1385) ANALYSIS 


Reaction 

Att 

A'y 

A'yTT 

A(1405) 

^ EV" 

0.083F0.004 

0.0007±0.0004 

0.658T0.012 

A(1405) 

^ E+tt" 

0.088±0.005 

0.0038F0.0009 

0.013T0.002 

A(1405) 

^ Ay 

0.008 ±0.003 

0.946F0.028 

0.098T0.009 

A(1405) 

^E^y 

0.585F0.019 

0.380F0.015 

0.837T0.023 

E(1385) 

^ Att 

0.905F0.010 

O.OllFO.OOl 

0.086T0.003 

E(1385) 

^ E+tt- 

0.050F0.002 

0.0018F0.0005 

0.00564T0.0008 

E(1385) 

—> Ay 

0.012F0.002 

1.309F0.022 

0.105T0.006 

E(1385) 

^E°y 

0.548F0.016 

0.24T0.01 

0.99T0.02 

A(1520) 

^ Ay 


1.388F0.027 

0.0010T0.0007 

A(1520) 

^E^y 


0.087F0.006 

0.586T0.016 

A(1520) 



0 

0.0099T0.0016 

A(1520) 

^E%° 


0.0006T0.0004 

0.681T0.014 


TABLE III: Acceptances (in units of 10“^) for the channels 
used in the calculation of the branching ratios. Here x%igh = 
y^ow ~ 3.841. The uncertainties are statistical only. 


Wigner of the form 


flMl 2MMoT{q) 

^ (M 2 - + M 2 r 2 (g) ’ 

Mq f Xq + Qq 


r(9) = ro ^ 


<70 


M V^o+9^ 


<7 = 


V'(M2 -Ml- M2)2 - 4M2M2 


2M 


go = 


V(M2 -Ml- M2)2 - 4M2M2 


2Mn 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


where Mq is the peak position of the resonance, Xo=0.35 
GeV and Tq is the width. For the E* ^ Att^ transi¬ 
tion, 1 = 1. We tried both first order and second order 
polynomial background parameterizations. The system¬ 
atic uncertainty in the yield extraction due to the choice 
of background function was about ±1%. The mass and 
width of the E(1385) were found to be 1.3860 GeV and 
0.03988 GeV. For the Ay channel (Fig. IT^). we used 
two relativistic Breit-Wigners (one for the E(1385) and 
one for the A(1520)) plus a polynomial background. The 
masses and widths were fixed to be those found from the 
fits to Figures [ni4 andll4L). 

From FigUV it is clear that we were not able to re¬ 
solve the A(1405) and the E°(1385), therefore in or¬ 
der to find the number of A(1405)’s (tia) we look at 
the events for which neither the y nor the 7r° hypoth¬ 
esis is satisfied (Fig. I12L>). This isolated predominantly 
A(1405) ^ E°7r° events, since the E(1385) ^ S°7r° decay 
is forbidden by isospin. We parameterized the A(1405) 
line shape using the M-matrix formalism for S-wave 
scattering below the KN threshold. The M-matrix is re¬ 
lated to the S-wave transition matrix T according to 


To obtain the yields we fitted the hyperon mass dis¬ 
tributions between 1.25 GeV and 1.75 GeV. The yield of 
Att^ events is extracted by fitting the data in Fig. If 4I \ 
with a polynomial background and a relativistic Breit- 


r = Q5(M-7Q)-igv (14) 

where Q is a diagonal matrix containing the relative E^tt'^ 
momentum q and KN momentum k |29j. Note that be- 
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Reaction 

Yield 

Estimated counts 

373.8±34.0 

Raw counts 

2878.3±77.4 

A(1405) ^ E°y 

0.45±0.17 

A(1405) ^ Ett 

95.7±9.5 

E(1385) ^ Ett 

10.4±1.0 

E(1385) ^ Ay 

0.87±0.21 

Corrected tt'^ counts 

2770.9±78.0 

Raw y counts 

100.2±15.4 

E(1385) ^ Att'’ 

35.0±1.0 

E(1385) ^ E+tt” 

0.38±0.3 

A(1405) ^ A°y 

0.85±0.27 

A(1405) ^ E°y 

0.29±0.11 

A(1405) ^ Ett 

2.47±0.25 

Corrected y counts 

61.2±15.4 


TABLE IV: Breakdown of statistics for the A 7 and Att® chan¬ 
nels. The errors are statistical only. 


low the KN threshold, the latter is purely imaginary. 
The matrix M is expanded relative to the KN threshold 
Et = Mm according to 


M{E) = M{Et) + \r{Q\E) - Q\Et)) 


Mil 

_ Mi2 


Mi2 
M 22 ’ 

k O' 

0 q ’ 


R = 


Rrn 

0 Rstt 


(15) 

(16) 
(17) 


result still needs corrections for S'*"?!" contamination 
and the contribution to the numerator from the reac¬ 
tion A(1405) ^ Ay. Based on the measured 27±8 keV 
radiative width^^, we assume that the leakage of the 
Ey channel into the 7 region is small relative to the Ay 
signal and that the leakage into the 7 r° region is small 
compared to the Att^ signal. The formula for the accep¬ 
tance corrected branching ratio is 


R 




An^ 


1 


An^A^(Ay) - An-),AJ(Ay) 


X 




jdTiTt 

2 



-An^ 



dStt 

^Att 

2 



(19) 


- N^{A* E+tt-) - N^{A* E°7r°) 

-A,(A*^E°y)-A^(A*^ Ay), (20) 

rij - N^{A* E+tt-) - N^{A* E°7r°) 

-A-,(A*^EOy)-A-,(A*^Ay), (21) 


where n-y (n,r) is the measured number of photon (pion) 
candidates and the remaining terms are correc¬ 

tions due to leakage from the A(1405). The acceptance 
for the individual pion (photon) channels are denoted 
as A?(E+ TT’ -), TT )) and so on. For example, 

A^(A 7 r) denotes the relative leakage of the Att channel 
into the Ay channel. Table ITTTI lists the values of these 
“acceptances”. 

The corrections depend on an estimate of the number 
riA of A(1405)’s in the data set. They are 


The amplitude for elastic scattering in the E° 7 r° channel 
is given by 

T22 = _ qiMii + \k\) _ 

(Mil + \k\){M22 - tq) - Mf^- ^ ’ 

Below Et, the Ett mass spectrum is proportional to 
|722p/g- Fig. inp shows the M-matrix parameteriza¬ 
tion fit to hyperon mass spectrum. A relativistic Breit- 
Wigner form is included to account for the leakage of 
the E(1385) —> Att^ channel into the high missing mass 
squared region. A second relativistic Breit-Wigner is 
used for the A(1520) contribution. The mass and width 
of the A(1520) were found to be 1.520 GeV and 0.022 
GeV. We used a second-order polynomial for the remain¬ 
ing background beneath the peaks. The matrix elements 
at threshold and the effective ranges were determined 
from the fit to be Mii(Et) = 1.314, Mi 2 (Et) = —1.063, 
M 22 {Et) = 0.686, R-rm ~ 9.543, and = —28.89. We 
find 328±36 A(1405) counts and 245±37 E(1385) counts 
in the hyperon mass region 1.34-1.43 GeV. The reduced 
for the fit was 0 . 866 . 

Although the 7r° leakage into the y channel is the 
dominant correction to the branching ratio, the final 


N^(A* ^ Ay) 
N^{A* E°y) 
N^{A* E°7r°) 

N^(A* E+tt") 


A^(A^)R(A* Ay)nA 

AA,(E07r0) + AA,(E+^-) 

A^{YP-i)R(A* ^ E0y)nA 
AA,(E07r0) + AA,(E+^-) 
AA(E07r0)nA 

AA,(E07r0) + AA,(E+^-) 
Ai^(E+7r“)nA 
AA,(E07r0) + AA,(E+^-) 


,( 22 ) 

,(23) 

,(24) 

,(25) 


and similarly for the pion channel. Here isospin sym¬ 
metry is assumed such that R{YPtt^) = i?(E+ 7 r“) = 
i?(E“ 7 r+) « 1/3 for the A(1405) decay channels. The 
subscript “y7r” refers to those events for which both a 
pion and a photon are missing or those events leaking into 
the “y7r” region due to the tail of the 7r° peak (this is why 
the E''' 7 r“ contamination must be included in the denom¬ 
inator, although the leakage for this channel is small). 
Table HVI lists the yields for the various channels of the 
E(1385) decays. The hyperon mass range was 1.34-1.43 
GeV. The reaction yp ^ AK*'^ causes a smooth back¬ 
ground underneath the E(1385) peak in FiguresEK and 
EP that is well parameterized by the second order poly- 
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nomial fit. Hence it has not been explicitly included Ta¬ 
ble |IV| The largest background in the 7 channel is due 
to leakage of the 7 r° tail into the 7 missing mass squared 
region. 

After subtracting the background contributions enu¬ 
merated in table nvi there were 61.2 ± 15.4 counts con¬ 
sistent with E°(1385) ^ A 7 and 2770.9 ± 78.0 counts 
consistent with E°(1385) ^ Att^. After correcting for 
the relative acceptance of the two channels, we obtained 
a branching ratio, of 


r(E°(1385) ^ A 7 ) 
r(E0(1385) ^ AttO) 


1.53±0.39(stat)%. 


(26) 


The branching ratio result for the E(1385) depends on 
how well we understand the tail of the 7 r° peak near the 7 
peak. Fig. ITTt l and CTp) shows the comparison between 
the data and the Monte Carlo for the reaction Y* AX 
for the 1.34-1.43 GeV hyperon mass region. The excess 
of counts above the 7 r° peak correspond to E ^ E° 7 r°, 
where E° —> A 7 . Although the Y A 7 decay is not com¬ 
pletely separated, a clear enhancement near zero missing 
mass squared can be seen above the tail clearly indi¬ 
cating the presence of radiative events. The Monte Carlo 
predicts that the leakage accounts for about 30% of the 
raw photon yield in the |M^| < 0.01 GeV^ region. In 
order to assess the quality of the Monte Carlo in the 
tail, we looked at A(1405) —> E+tt” events for which the 
E+ subsequently decayed to p7r°. We chose this chan¬ 
nel because there are no channels that can distort the 
spectrum above the 7 r° peak, the E"*" radiative channel is 
rare (BR = 1.25 x 10“^), and has similar kinematics to 
the Att^ decay. We required the p7r~ invariant mass to be 
greater than 1.13 GeV (to eliminate the A(1116) from the 
sample). To identify the E+(1189) we require the pX in¬ 
variant mass (or, equivalently, the missing mass recoiling 
off the AT+tt” system) to be in the range 1.17-1.206 GeV. 
We performed kinematic fits on these events with vertex 
and four-momentum conservation constraints. Explicitly, 
the constraint equations are 


/ = 


Bbeam T Br Bit By;+ 

Pbeam Ptt Pp Ps+ 

iv - y^)pl - {z - z^)py 
{x - x^)pI - (z- z^)p% 

{y - yp)PK -{z- zp)p\ 

{x - Xp)pl^ -{z- Zp)p^ 


= 0 . (27) 
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FIG. 17: Comparison between data and Monte Carlo results 
for the reactions and Y* AX (top histograms) Y* —> E'^ 7 r“ 
(bottom histograms) after kinematic fitting has been per¬ 
formed. The points with error bars are the data and the 
curves are the MC results. Histrograms B and D have the ver¬ 
tical scales expanded by a factor of ten. In B the solid curve 
on the left is the E* —> A 7 simulation, the central dashed 
curve is the E* ^ Att'^ simulation, the isolid curve on the 
right is the A(1405) ^ simulation. In A the curve is 

the sum of the three. In C and D the E+tt” data and the 
E(1385) —> Att® Monte Carlo distribution have been scaled 
to agree with the peak height of the 7 r° in the Y* —> AX 
distribution from the data set. 


-bO.Ol GeV'^ region is about ^ 19%. Scaling the leakage 
of the Att^ channel into the 7 region by a factor of 1.19 
reduces the branching ratio from 1.53% to 1.36% for a 
relative change of about —11%. More importantly, com¬ 
paring with UZp shows a clear enhancement at zero 
missing mass present for the latter case not in evidence 
for the former case. The negative systematic error will 
be increased by 11 % in quadrature. 


A(1520) ANALYSIS 


The missing mass squared distribution for the reaction 
chain Y* —> E+ 7 r“, E+ ^ p{X) is shown in Fig. 117C 
and CTp for hyperon masses in the 1.38-1.45 GeV mass 
region. We used the four-vector for the E+ obtained 
from the fit, with less than 0.5% probability of exceed¬ 
ing The Monte Carlo result (dashed histograms) for 
the A(1405) ^ E+tt" reaction agrees very well with the 
data down to about zero missing mass squared. The dis¬ 
crepancy between the MC and the data in the —0.01 - 


For the A(1520) analysis we calculated the radiative 
branching ratio relative to the E° 7 r° and the E+ 7 r“ chan¬ 
nels. The hyperon mass cut used to identify the A(1520) 
was 1.49-1.55 GeV. From the fit to the histogram shown 
in Fig. El we obtained = 32.5 ± 8.2. To identify the 
E^tt*' channel we used events for which neither the 7 nor 
the 7 r° hypothesis is satisfied. The ground-state A is a 
decay product in the E° 7 r° (14%) and Atttt (10%) chan- 
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FIG. 18: Sample fit of the AX mass distribution for missing 
mass squared in the 0.018-0.075 GeV^ range. 


Reaction 

Yield 

A(1520) ^ E+tt" 

5290±124 

A(1520) ^ EV*" 

202.8±16.7 

A(1520) ^ A 7 

0.05±0.01 

Gorrected counts 202.8±16.7 

Raw 7 counts 

32.5±8.2 

A(1520) ^ E^tt^ 

0.09±0.01 

Gorrected 7 counts 

32.4±8.2 


TABLE V: Breakdown of statistics for the A(1520) analysis. 
The errors are statistical only. 
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FIG. 19: Isolation of Y* —>■ events. A) and B) show 

the K^-k~ vertex distributions. C) is the missing mass for the 
reaction qp —> K^ti~X. D) is the hyperon mass distribution 
for events satisfying the E'^ identihcation cut (see text). 


nels. In order to simplify the calculation for the branch¬ 
ing ratio we require the missing mass squared to be in the 
range between m^o and 0.075 GeV^ (« 4to^o, the two- 
pion threshold). This isolates the S° 7 r° channel. The hy¬ 
peron mass distribution in the A(1520) region with this 
additional cut applied is shown in Fig. The fit is a D- 
wave {I = 2) relativistic Breit-Wigner plus a polynomial 
background. We tried both first-order and second-order 
polynomials; the results for the yield differed by ± 1 . 6 %. 
The leakage of one channel into the other is neglible and 
applying the correction does not change the result. Due 
to the low acceptance for events containing A’s, the raw 
number of E° 7 r° counts is only a factor of 6 larger than 
the radiative signal and the technique relies on isolat¬ 
ing a channel for which two particles (7 and 7 r°) are not 
detected. We also looked at A(1520) ^ E+tt” events 
for which the acceptance is higher. The same particle 
identification and vertex cuts used for the previous anal¬ 
ysis were applied with some modifications. We required 
that the p'K~ invariant mass be greater than 1.13 GeV 
to cut A(1116) contamination. The primary vertex was 
determined using the K'^ and tt~ tracks. The z-position 
and X- and y-positions for these vertexes are shown in 
Fig. 1191 4 andll9B. respectively. A prominent E+(1189) 
peak shows up in the missing mass recoiling against the 
K'^ and the TT (Fig. n^ti). The hyperon mass spec¬ 
trum for those events in the range 1.165-1.215 GeV about 
the E"*" peak are shown in Fig. cnp. The curve is a fit 


to the A(1520) region using a D-wave relativistic Breit- 
Wigner with a second-order polynomial background. In 
the region between 1.49 GeV and 1.55 GeV we obtain 
riA* = 5290±124 (the acceptance of CLAS is much larger 
for this channel than the others due to the larger 7 r“ mo¬ 
mentum) . The yields for these two reactions are listed in 
table m As can be seen from the numbers in the table 
the leakage of each channel into the other is negligible. 
None of the generated A(1520) ^ E+ 7 r“ events satisfied 
the selection criteria. There is no Att® leakage since this 
channel is forbidden by isospin. 

We obtained a raw branching ratio of rijfn-yn = 16.0 ± 
4.3%. Gorrecting for acceptance, the branching ratio is 


r(A7) 

r(E07r0) 


A^(A7) 'IT'-y-K 


7.9 ±2.1%. 


(28) 


The acceptances used in this calculation are listed in 
table EH To obtain the branching ratio r(A 7 )/F 7 ’ 07 ’ 
we scale this result by the branching fraction of 14% for 
the E° 7 r° channel (assuming isospin symmetry) to obtain 
1.10±0.29%. The acceptance for the A(1520) ^ E+ 7 r“ 
channel was 1.66 ± 0.06%. We obtain 1.01 ± 0.26% for 
the radiative branching ratio. The results for the two 
channels agree after acceptance corrections. 

If contamination due to the A(1520) ^ E °7 channel is 
present the branching ratio for the A 7 channel acquires 
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X^ow Xhigh R(%) S(1385) R(%) A(1520) 


2.706 

2.706 

1.68±0.41 

1.20T0.29 

3.841 

3.841 

1.53±0.39 

1.10±0.29 

6.635 

6.635 

1.58±0.40 

1.13±0.33 

2.706 

6.635 

1.38±0.38 

1.25±0.33 

3.841 

6.635 

1.36±0.30 

1.06T0.34 


TABLE VI: Dependence of the S(1385) ^ A 7 and A(1520) —> 
A 7 branching ratio on the choice of Xhigh,low cuts. 


largest theoretical estimate for the E °7 channel relative 
to the A 7 channel of 0.153 from R. Bijker, F. lachello, 
and A. Leviatan^^, the partial width is reduced to 478 
keV, which is an insignificant change. 

For the A(1520) decay, we obtained a branching ratio 
of 1.10 ± 0.29(stat)to;o4(sys)9o using the E° 7 r° channel 
and 1.01 ±0.27% using the S+ 7 r“ channel. The weighted 
average gives a branching ratio of 


a small correction term: 


R{Ax) 


A (S°7r°) 

±R(E°7) ( -^7^ ~ \^9) 


where i?(E° 7 ) is the branching ratio to the E °7 channel 
and i?(E° 7 r°) is the branching ratio to the E° 7 r° channel. 
Using the largest theoretical estimate for the E °7 radia¬ 
tive width of 293 keV from Warns, et aZ.j^, we obtain a 
correction of ±0.01%. Therefore this contamination can 
be neglected. 


RESULTS 

To check the sensitivity to the conhdence limits used, 
i?^^o was calculated with 1%, 5% and 10% probability 
for accepting a channel and 99%, 95% and 90% prob¬ 
ability for rejecting a channel. Table ED lists the cor¬ 
rected branching ratios as a function of the Xhigh low 
cuts. The third column in Table ED gives the 2(1385) 
results. The results were very stable, varing from ±0.15 
(10%, 90%) to —0.17 (5%, 99%). These values were used 
as estimates of the systematic errors. The value for the 
branching ratio is 1.53±0.39(stat)lo :[ 7 (sys)%, where the 
second uncertainty reflects the variation in the branching 
ratio as a function of the choice of cuts. 

We add the 11% relative error (i.e. —0.17% absolute) 
that could result from underestimating the tail of the 
TT^ response to the negative systematic error and quote 
a branching ratio of 1.53 ± 0.39(stat)to;24(sys)%. The 
positive systematic error reflects the range of values we 
obtained for the various estimates for the branching ratio. 
If we neglect the small (unmeasured) contribution due to 
the E °7 channel, the S°(1385) ^ A 7 partial width is 
given by 

= 479±120(stat)l^Jo(sys) keV, 

(30) 

using Ftot = 36 ± 5 MeV and = 0.135 ± 0.011, 
the branching ratio of the Stt channels relative to the 
Att® channel!^. The errors on Ftot and R ^^ are in¬ 
cluded in the systematic error for r(A 7 ). If we use the 


r(A(1520) ^ Ay) ^ 1 07 ± 0.29(stat)t° il(sys)%. (31) 
i tot 

Table IVll lists the branching ratios for various combina¬ 
tions of kinematic fitting cuts. There is no obvious 
dependence on the choice of cuts. To determine the sys¬ 
tematic error in the measurement using the channel 
to normalize, we used the range of branching ratio val¬ 
ues obtained for different choices of cuts. Using a full 
width of 15.6±1 MeV|i3, we obtain a partial width of 
167±43(stat) 1 * 112 ( 85 ^ 8 ) keV. The error on the full width is 
included in the systematic error for F(A 7 ). The A(1520) 
result is compatible with the Mast et al. result and the 
Antipov et al. result [l^ but disagrees with the Bertini et 
al. result 0. Together, our result and those of Mast et 
al. and Antipov et al. exclude the bag models listed in 
Table ID 

The E°(1385) ^ Ay channel has never been measured 
before. The result is roughly 2-3 times larger than all 
of the existing model predictions except for HByPt[i^. 
Table^reveals that the model predictions for the A ^ py 
transition are also about 50% low. Sato and Leej3 ij 
showed that much of that discrepancy could be accounted 
for by the inclusion of non-resonant meson-exchange ef¬ 
fects. They find a width of 530±45 keV, about 80% of the 
experimental value. Lu et al. reproduced the A ^ py 
data using a chiral bag model calculation with a rela¬ 
tively small bag radius of 0.7 fm. About 40% of the 
transition was due to the pion cloud. These calculations 
suggest that mesonic effects could account for the dis¬ 
crepancy between the model predictions and our result 
for the S°(1385) radiative transition. 
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